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We present a method of performing high-speed rotational anisotropy nonlinear optical harmonic 
generation experiments at rotational frequencies of several hertz by projecting the harmonic light 
reflected at different angles from a sample onto a stationary position sensitive detector. The high 
rotational speed of the technique, 10^ to 10^ times larger than existing methods, permits precise 
measurements of the crystallographic and electronic symmetries of samples by averaging over low 
frequency laser power, beam pointing, and pulse width fluctuations. We demonstrate the sensitivity 
of our technique by resolving the bulk four-fold rotational symmetry of GaAs about its [001] axis 
using second harmonic generation. 


Rotational anisotropy nonlinear harmonic generation 
(RA-NHG) is an all-optical technique for determining 
the crystallographic and electronic symmetries of solids 
[1]. It has been successfully applied to study a variety 
of phenomena, including lattice reconstruction on semi¬ 
conductor thin film surfaces mg and metal/electrolyte 
interfaces mil], magnetic ordering on metal surfaces and 
buried multilayer interfaces IHHIS], magnetic ordering in 
bulk semiconductors m and multiferroics [EHIH], and 
molecular self-assembly on dielectric substrates [19]. In 
the RA-NHG technique, a laser beam is focused onto a 
sample and the intensity of light reflected or transmitted 
at a higher harmonic of the incident frequency is mea¬ 
sured as a function of the angle (f between the scattering 
plane and some crystalline axis parallel to the surface of 
the sample while a fixed angle of incidence is maintained 
[inset Fig. [^. By selectively measuring different com¬ 
binations of incoming and outgoing light polarizations, 
typically linearly polarized either in (P) or out of (S) 
the scattering plane, the independent non-zero elements 
of the nonlinear optical susceptibility tensor of the sample 
can be deduced. By Neumann’s principle, this serves to 
identify the crystallographic and electronic point group 
symmetries of the sample under study [Tj. 

RA-NHG experiments are typically carried out by ro¬ 
tating the sample while keeping the scattering plane 
fixed. An alternative approach in which the scattering 
plane is rotated while the sample is kept fixed has also 
recently been demonstrated, which is advantageous when 
studying small samples and when low temperatures or 
externally applied fields (e.g. magnetic, electric, strain) 
are desired [201 ED- In both methods, (f> is stepped and 
held at discrete values for a finite wait time while the high 
harmonic intensity is recorded. A major drawback of this 
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method of data acquisition is the long time required (of 
order 10^ sec) to complete a full 360° sweep of (f. This 
introduces a large amount of low frequency noise in the 
data from temporal fluctuations in laser power, beam 
pointing, and pulse width, limiting the ultimate preci¬ 
sion with which point group symmetries can be deter¬ 
mined. It also precludes accurately measuring temporal 
changes in the RA-NHG response of a sample that occur 
on the timescale of a 0 rotation period, which may arise 
from processes such as charge migration or impurity ad¬ 
sorption on the sample surface under vacuum conditions 
[22l 1^ . In this Letter, we demonstrate a novel high¬ 
speed RA-NHG setup using a position sensitive detector 
that rotates 0 at a frequency 10^ to 10^ times higher 
than conventional methods [20|, which can in principle 
be increased to even higher frequencies. 

A schematic layout of the high-speed RA-NHG setup 
is shown in Fig. Unless otherwise stated, all optical 
components are standard and were obtained from Thor- 
labs, Inc. The setup operates in the following way: A 
circularly polarized laser beam, produced from a linearly 
polarized beam using a quarter-wave plate (QWP), first 
passes through a linear polarizer (LPl) and a fused sil¬ 
ica binary phase mask (PM; Tessera Technologies, Inc.) 
that are mechanically co-rotating about the beam axis. 
The polarizer produces a beam of constant intensity with 
a rotating linear polarization that is maintained either 
parallel or perpendicular to the PM groove direction de¬ 
pending on whether an Sin or Pin polarization geometry 
at the sample is required. The phase mask splits the 
beam into many diffracted orders, each retaining the po¬ 
larization set by LPI and sweeping out a cone as the 
PM rotates. For our measurements, only the +I order 
is utilized while all other orders are blocked. Next, a 
collimating lens (LI) and an achromatic objective lens 
(L2) form a telescope that focuses the beam onto the 
surface of the sample, thus creating a rotating scatter¬ 
ing plane. The fundamental and higher harmonic beams 
reflected from the sample return via the diametrically 
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FIG. 1: Schematic layout of the high-speed RA-NHG setup. Laser light is incident from the left and passes through a series 
of optics labeled as follows: quarter-wave plate (QWP), linear polarizers (LP1,2), binary phase mask (PM), collimating lens 
(LI), dichroic mirrors (DMl-3), achromatic objective lens (L2), spectral filter (SF), and GGD camera (GGD). LP2 is mounted 
on one side of a rotating wheel while the other side has a through-hole. Only the +1 diffracted order from the PM used in 
our experiment is shown. The polarization of the fundamental (red) and higher harmonic (blue) light is illustrated at various 
points along the beam path. The PM and two LPs are mounted on a motorized rotation stage (not shown) and are rotated 
co-axially (solid curved arrows). The dashed arrows show the path traced by the beam on LI, L2, DMl-3, SF, and GGD during 
the rotation. Top left inset shows the simplified geometry of a reflection-based RA-NHG measurement. 


opposite side of L2 upon which they are re-collimated. 
They then pass through a linear polarizer (LP2) that 
is co-rotating with both LPl and PM, which is used 
to select either the Sout oi* Pout polarization component 
of the reflected beam. The desired harmonic is finally 
isolated using tilted dichroic mirrors (DM2,3; Semrock, 
Inc. BrightLine series) and spectral filters (SF) that di¬ 
rect it onto a cooled electron-multiplying CCD camera 
(CCD; Andor Technology Ltd.). As LPl, PM, and LP2 
rotate, the reflected beam sweeps out a circle on the two- 
dimensional active area of the CCD, thereby recording 
the angular ^-dependence of the NHG signal as a posi¬ 
tion dependent intensity on the CCD. To eliminate any 
possibility of mechanical phase slip between LPl, PM, 
and LP2 during rotation, they are coupled to a common 
axle driven by a single 12 V DC electric motor with a 
pulse-width modulating power controller. The need to 
rotate only three lightweight optical elements allows for 
rotation at high speeds (^ 4 Hz or more), which is greater 
than previous setups that rotate the detector [20l |2T] by 
a factor of 10^ to 10^ and much faster than the character¬ 
istic timescale for laser fluctuations and drift. Moreover, 
because the value of 0 read from the CCD image is rigidly 
tied to the angle of the PM, the experiment is insensitive 
to random phase slips of the motor. 

The key component of our high-speed RA-NHG setup, 
enabling the use of a stationary detector, is the triple 
dichroic mirror periscope (DMl-3). This consists of three 
identical long-pass dichroic mirrors with a cut-on edge 
wavelength between that of the fundamental and desired 


harmonic, each mounted at 45° with respect to the di¬ 
rection of beam propagation, as shown in Fig. To 
illustrate the function of this periscope, we first study its 
effect on the fundamental beam by measuring the inten¬ 
sity transmitted through DM1, DM2, and both mirrors 
as a function of the incident linear polarization angle cj) 
[Fig.m: a)]. As (/) changes, the polarization incident on a 
dichroic mirror that is rotated by 45° about the vertical 
axis varies continuously between S and P with a 180° 
periodicity. The transmittance through DM1 therefore 
varies sinusoidally in accordance with the Fresnel equa¬ 
tions, as plotted in Fig.j^b). In order to compensate for 
this unwanted effect, the beam is sent through a second 
dichroic mirror (DM2) that is tilted by 45° about the 
horizontal axis, which has a transmittance curve that 
is 180° out of phase with the first mirror [Fig. ic)|. 
The combined effect of DM1 and DM2 is to completely 
eliminate any net polarization dependence of the trans¬ 
mittance [Fig. I^d)] and thereby maintain a constant <p- 
independent incident intensity on the sample. The com¬ 
bination of DM2 and DM3 acts in a similar way in re¬ 
flection geometry to eliminate any net polarization de¬ 
pendence of the outgoing high harmonic beam. We note 
that the acquisition time for each curve shown in Fig.j^is 
approximately 10 minutes. On this timescale, character¬ 
istic low frequency laser fluctuations are clearly observed 
through the deviations of the data away from the fitted 
curves and the mismatch of the 0 = 0° and 360° data 
points. It is precisely these types of fluctuations that are 
averaged out by the high-speed RA-NHG setup. 








3 



0 40 80 120 160 200 240 280 320 360 


Incident Polarization Angle (|) (deg) 

FIG. 2: Demonstration of polarization compensation by a 
dichroic mirror pair, (a) Polarization-dependent transmit¬ 
tance curves of DM1 and DM2 were measured using a photo¬ 
diode positioned after DM2 with both PM and LI removed. 
Data were acquired by stepping the polarization angle 0 in 6° 
increments and averaging for ~ 10 sec per data point. The 
transmitted intensity of 800 nm light through (b) only DM1, 
(c) only DM2, and (d) both DM1 and DM2 are displayed as 
a function of cj). Also shown are fits to cos cj) for the individual 
mirrors and a constant for the mirror pair. Fluctuations of 
the data about the fitting functions are due primarily to low 
frequency laser power fluctuations during acquisition. The 
error bars show the RMS deviation of the data from the fits. 


To characterize the performance of our RA-NHG 
setup, we measured the rotational anisotropy of second 
harmonic generation (SHG) from the (001) crystal sur¬ 
face of GaAs. For the incident light source we used a 
regeneratively amplified Ti:sapphire laser system (KM- 
Labs Wyvern 1000) that produces 60 fs optical pulses 
(broadened to ^ 150 fs at the sample by dispersion) with 
800 nm center wavelength and operates at a 10 kHz rep¬ 
etition rate. The dichroic mirrors were selected to have a 
cut-on wavelength of 635 nm. At a rotation frequency of 
{l/360°)d(l)/dt 4 Hz, the angular spacing between suc¬ 

cessive pulses is ^ 0.14°, which is much smaller than the 
angular subtense of each beam spot on the CCD (^ 9°) 
and therefore allows the incident beam to be treated as 
quasi-CW. Raw CCD images of the RA-SHG data (400 
nm wavelength) from (001) GaAs obtained in a Pin-Sout 


geometry are shown in Fig.j^a) for a series of short con¬ 
tinuous camera exposure times. A ring of intensity with 
alternating nodes and antinodes separated by 45° is al¬ 
ready clearly visible after 1 sec, which demonstrates a 
rotationally anisotropic SHG response of the sample. 

The CCD images can be cast into a form more con¬ 
ducive to quantitative analysis by performing a radial 
integration (after background subtraction) with respect 
to an origin centered on the ring of intensity and binning 
the data into fixed <p azimuthal intervals. Figure |^b) 
shows the integrated intensity of both the 5 sec and 20 sec 
images as a function of 0 on a polar plot, which exhibit 
four lobes of equal amplitude consistent with previously 
reported RA-SHG measurements on the (001) GaAs sur¬ 
face [2]. Quantitatively, the data are also well described 
by the expression |Acos(20)p expected under Pin-Sout 
geometry, where A is a linear combination of bulk and 
surface electric dipole induced SHG susceptibility tensor 
elements derived from the C 2 v symmetry group of GaAs 
[2] . A closer examination of the data reveals a high degree 
of overlap between data points taken at different expo¬ 
sure times, which demonstrates the high level of precision 
achievable using our technique and shows that the minor 
deviations of the data away from the ideal |Acos(20)p 
form are due to systematic sources of error such as mis¬ 
alignment and sample surface quality rather than statis¬ 
tical sources of error such as laser fluctuations. 

In summary, we have demonstrated a novel technique 
based on position sensitive detection to perform RA- 
NHG experiments at rotational frequencies 10^ to 10^ 
times higher than existing methods. The frequency can 
in principle be pushed much higher using a more pow¬ 
erful motor m, which would further suppress low fre¬ 
quency noise arising from laser fluctuations and allow 
increasingly sensitive measurements of point group sym¬ 
metries. This technique can be applied in both reflection 
and transmission geometries and can be extended over 
a wide photon wavelength range, limited primarily by 
the spectral range of the CCD. It is also amenable to 
ultrafast time-resolved pump-probe RA-NHG measure¬ 
ments through the use of short laser pulses and the in¬ 
troduction of a pump beam, which can be used for ex¬ 
ample to separate electronic relaxation pathways |25] or 
to study photo-induced symmetry breaking phase tran¬ 
sitions [26] . The time resolution of such experiments can 
be improved by replacing the transmissive objective L2 
with a reflective Cassegrain objective to reduce dispersive 
pulse broadening, and sensitivity to small pump-induced 
changes in NHG intensity can be enhanced with a lock-in 
camera. 
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FIG. 3: RA-SHG data collected from (001) GaAs using the high speed RA-NHG setup, obtained with 800 nm incident light 
under a Pin-Sout polarization geometry and a fluence of ~ 5 mJ/cm^. (a) GGD images with 1, 5, 10, and 20 sec exposure times. 
The GaAs crystal was aligned such that 0 = 0° corresponds to the [100] axis, (b) Polar plot of the RA-SHG intensity obtained 
by performing a radial integration of the 5 sec and 20 sec data in panel (a) over the interval bounded by the concentric white 
circles, with an azimuthal bin size of A0 = 6°. A fit to the expected functional form derived in Ref. [2] is overlaid. We note 
that because the camera exposure times are not locked to integer multiples of the 0 revolution period, there will be a partial 
range of angles that acquire one extra revolution’s worth of signal. This effect can be corrected by reducing the overshot data 
points by a factor of 1/iV where N = (exposure time x rotation frequency) is the total number of revolutions acquired. Such 
corrections are only important for short exposure times and were applied to the 5 sec data between 98° and 134°. 
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